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We clarify the relation between the improvement of dispersion relations in the fermion sector of 
00 ■ lattice regularized QCD and the improvement of bulk thermodynamic observables. We show that 

' in the infinite temperature limit the cut-off dependence in dispersion relations can be eliminated 

' up to 0{a") corrections, if the quark propagator is chosen to be rotationally invariant up to this 

order. In bulk thermodynamic observables this eliminates cut-off effects up to the same order at 
(-H ' vanishing as well as non-vanishing chemical potential. We furthermore show, that in the infinite 

' temperature, ideal gas limit the dependence of finite cut-off corrections on the chemical potential 

is given by Bernoulli polynomials which are universal as they do not depend on a particular dis- 
cretization scheme. We explicitly calculate leading and next-to-leading order cut-off corrections for 
' some staggered and Wilson fermion type actions and compare these with exact evaluations of the 

free fermion partition functions. This also includes the chirally invariant overlap and domain wall 
fermion formulations. 
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r-| ; I. INTRODUCTION 

Numerical studies of lattice regularized QCD face the problem of lattice discretization errors which complicate the 
^ ' extraction of physical results in the continuum limit. This problem has been addressed ever since the formulation 
\ of a systematic improvement scheme for gauge theories by Symanzik [1]. Aside from the now widely used Symanzik 
OO . improved gauge actions it also led to tree-level improved actions for the fermion sector of QCD like the Naik [2] and 
p4 [3] staggered fermion actions as well as the clover [4] and (truncated) perfect actions for Wilson fermions [5, 6]. 

The reduction of cut-off effects plays a particularly important role in studies of QCD thermodynamics where the 
relevant observables, like energy density or pressure, are dimension four operators and the numerical signal for these 
^5 observables thus drops like the lattice spacing to the fourth power. This forces one to perform calculations on rather 
coarse lattices on which discretization errors can be significant. In fact, this has been observed early on in studies 
of the bulk thermodynamics of SU(3) gauge theories [7] and the efficiency of improved actions for studies of gauge 
theories at high temperature has been demonstrated [8]. 

In calculations on lattices with finite temporal extent N^, as required in thermodynamic studies, cut-off effects 
obscure strongly the approach to the infinite temperature, ideal gas limit. Improvement schemes for thermodynamic 
' calculations therefore have been developed and tested quite successfully in this limit for both the gauge as well as the 
fermion sector of QCD. The qualitative features of the cut-off dependence of thermodynamic observables present in 
this limiting case have been found to carry over to numerical calculations at finite temperatures and are significant 
even at temperatures as low as a few times the transition temperature to the quark-gluon plasma phase. 

In this paper we will concentrate on a discussion of cut-off effects in the fermionic sector of bulk thermodynamic 
quantities that are given by derivatives of the QCD partition function with respect to temperature, T, or quark 
chemical potential, /i. In particular the latter aspect has so far not been analyzed systematically. Some results for the 
/z-dependence of cut-off effects have been obtained for the Naik and p4 staggered fermion actions [9, 10], truncated 
perfect actions [11] and for overlap fermions [12]. In fact, in the latter case the question was raised whether the 
introduction of a non-zero chemical potential in the overlap formalism could spoil the leading order cut-off dependence 
of thermodynamic quantities and could reintroduce divergences of the type discussed in early studies of QCD at finite 
density [13]. We proof here that this is not the case and give explicit results for the /i-dependence of cut-off effects 
in leading and next-to-leading order in a large- AT^ expansion of bulk thermodynamics. We show that a non-vanishing 
chemical potential only modifies the expansion coefficients but does not change the structure of the expansion, i.e. 
actions that are improved to 0{a^) at /i = remain improved to that order also for yu > 0. We will clarify the relation 
between the improvement of dispersion relations and the improvement of bulk thermodynamic observables and give 
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explicit results for the large-A',- expansion of the pressure calculated with staggered and Wilson type actions. Wc 
furthermore show that these considerations carry over to chirally invariant fermion formulations such as overlap and 
domain wall fcrmions. 

We will analyze the relation between rotational invariance of quark propagators, improved dispersion relations and 
bulk thermodynamic observables in the next two sections. We further exemplify these results for naive staggered 
fermions as well as two popular improved staggered fermion actions, the Naik [2] and the p4 action [3] in section IV. 
In section V we consider fermion discretizations of the Wilson type with general couplings on a hypercube of size 
(2a)''. We finally comment on thermodynamics with overlap and domain wall fermions in section VI and give our 
conclusions in section VII. Appendices A through D contain some useful formulas to make the paper self-contained. 



II. DISPERSION RELATIONS AND THE FREE ENERGY 



We will start our discussion of cut-off effects for free fermions which corresponds to the asymptotic infinite tem- 
perature limit of QCD by considering the simplest (naive) lattice discretization scheme used for the fermion sector of 
QCD. This will already demonstrate the main features of our analysis. For simplicity we also will suppress color and 
flavor factors, which are just multiplicative in the infinite temperature limit. 

The partition function for free fermions at temperature T and for non- vanishing chemical potential /i, defined on a 
lattice of size N^Nr, is given by 

Z{V, T,fi,m)= III dVixdVx e-^^ , (1) 

where {N^a)^ = V and N^a = 1/T are the volume and inverse temperature of the system; a denotes the lattice 
spacing and Sp is the Euclidean action for free fermions of mass m. 

In the simplest discretization scheme for fermionic actions derivatives are replaced by nearest neighbor differences 
on a four dimensional lattice, 

X \k=l J 

Here we also introduced the chemical potential through the usual exponential form [13]. The; Grassmann- valued 
fermion fields obey anti-periodic boundary conditions in the time direction. We note that the transformation tpx,x4 
e~^^'^tpx,x4, V'x,x4 — * e'*^*Vx,x4 leaves the path integral over the fermion fields invariant and shifts the entire fi- 
depcndence of the fermion action into the last time slice. The Euclidean action thus only depends on the chemical 
potential in units of the temperature, /x/T = ^aN^.. The action can be written in momentum space as ^ 

Sf = Y, ^{p)ilMP, I^Mp) (3) 
This leads immediately to the denominator of the fermion propagator 

4 3 

D{p,P4, m) = X] ^'^(^'^'4, fJ')D^{p,p4, /i) = ^ sin^(apfe) + sin^(ap4 - i^a) , (4) 
i/=i fe=i 

and allows to evaluate the partition function explicitly, 

z{v,T,^i,o) = llD^p,P4,^i) . (5) 
p 

Here the momenta take on discrete values, apk = 2'Knk/N„, rife = 0, ,±1, ...,±{N„/2 — 1), N„/2 and ap4 = 27r(n4 -|- 
l/2)/Ar,, n4 = 0, ,±l,...,±(Ar,/2-l), - Nr/2. 

For our general discussion of cut-off effects that arise from the finite lattice spacing given in units of the temperature, 
aT = 1/Nr, the spatial extent of the lattice is of less interest. We thus take the thermodynamic limit, A^^ — > oo. 



^ Here and in the following we consider massless fermions, with obvious generalization for the massive case. 
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which simpHfics the following considerations. We then obtain the pressure of a free fermion gas at temperature T and 
for chemical potential /k, 



d^av 1 

= 2/ TCTaT E ln[a;'(p)+sin2(27r(n + l/2)/Ar.-iH] , (7) 

with cj^(p) = sin^(apfe). Note that we reduced the integration and smumation intervals to half the Brillouin 

zone which eliminates a factor 16 arising from the 16 so-called doublers which emerge in the naive discretization. We 
thus normalize to a single quark flavor as we also will do for other discretization schemes in the following sections. In 
order to evaluate the sum over Matsubara modes appearing in Eq. (7) one uses the fact that this sum can be viewed 
as the sum over residues resulting from a contour integral over poles in the complex plane [14], 

J2fi^n)^ <f^J{r)hir,f,). (8) 

n 

The computational steps given in [14] can easily be generalized to non vanishing chemical potential by noting that 
the entire dependence is contained in the function /i(r, /x), 

M^>M) = ;p?rT^^> (9) 

the poles of which at 

rn = cx.p[ii2Trin+l/2)/Nr-ina)] (10) 
deliver the sum over the Matsubara frequencies via the relation 

sin^ {2TT{n + l/2)/Nr - i^^a)) = ~ (^rn- . (11) 
Carrying out the sum over the discrete set of momenta yields 

N-r/2-l 

J2 ln[sin2(27r(n + l/2)/7V^ - i^ia) + uj^{p}] = In (l + ze-^-^^^P)^ + In (l + ^-ie-^-«-E®) + const. (12) 

71=0 

where z = cxpipaNr) = cxp(/j,/T) is the fugacity and E(p} is the dispersion relation which is given by the pole of the 
fermion propagator, i.e. E{p) is obtained from the zeroes of Eq.(4) for jia = 0, 

D(p, E = ip4, 0) = <^ LO^{p) - sinh^iaE) = . (13) 

This emphasizes the particular role dispersion relations play in the study of bulk thermodynamics. 

Inserting Eq. (12) into Eq. (7) and subtracting the zero temperature part of the pressure, (Pa^)o = limAr^^oo Pa'^ 
eliminates ultra-violet divergences and allows to define the pressure in the same normalization as it is used in numerical 
calculations on the lattice, 

f_ ^ [Pa^-(Pa\]N^ 

= / d^ap [ln(l + ;sexp(-Af^a£;(p)))+ln(l + ^-^exp(-iV^a£;(p)))] . (14) 

t/— 7r/2 

Eq. (14) is the starting point for a systematic analysis of cut-off effects in bulk thermodynamic observables that 
are obtained from the logarithm of the partition hmction, P/T* = limv'^oo(^r^)~"'^ In Z, in terms of derivatives with 
respect to T or /x. An expansion of P/T** around 1/Nt = yields systematic corrections to the continuum ideal gas 
result which are given in terms of even powers of l/Nr- This arises after introducing yj = Pi/T = N^api as new 
integration variable in Eq. (14) and realizing that the dispersion relation is an even function in the momenta. We will 
discuss this expansion in more detail in the next section. 
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III. IMPROVEMENT OF DISPERSION RELATIONS AND BULK THERMODYNAMICS 
A. Generalized dispersion relations and the pressure 

We will outline here the arguments that lead to the observation that also for more general fermion actions than the 
one discussed in the previous section there exists a close relation between that branch of the dispersion relation that 
survives in the continuum limit and the large- TV,- expansion of P/T^. 

Similarly to the fermion propagator derived in the previous section for the naive discretization scheme of the 
fermion action one obtains results for general fermion actions of staggered or Wilson type that also may involve more 
complicated terms than the 1-link term used in the previous section. The propagator for free fermions of staggered 
type, D(p,p4), can generally be written as a polynomial^ in sin^(ap4), 

n 

D{p,P4) = "^i^P) sin''(aP4) , (15) 

where the coefficients di depend on the spatial momentum components. For Wilson type quarks the inverse propagator 
can instead be written as a polynomial in sin^(ap4/2), 

n 

D{p,P4) = Yl ^^(^ sin2^(ap4/2) . (16) 

i=0 

This is due to the Wilson term which introduces a cos(ap4) = 1 — 2 sin^(ap4/2). In both cases, rewriting D in terms 

of its n roots, (jJi{p), 

n 

D{p,k4)=dn{p)l[[sm\ak4)+ujf{p)] , (17) 

where k4 = p4 for staggered and k4 = P4/2 for Wilson quarks, immediately gives the n dispersion relations, 

smh^{aei{p)) = — sin^(afc4) = a'l(p) , i = l,..,n , (18) 

with Ci = Ei for staggered and = Ei/2 for Wilson type quarks. Although D is real, the roots can in general be 

complex. However, at least for small momenta the root which survives the continuum limit, loi to be definite, must 
be real in order to reproduce the continuum dispersion relation, E'^{p) — p^ . When we insert the right hand side of 
Eq. (17) in the partition function, Eq. (5), we find that each factor appearing in the product over the n distinct roots 
contributes in the evaluation of P/T^ as a term in a sum over the different branches, Ei, of the dispersion relation. 
Each term in this sum has the same structure as the integral appearing in Eq. (14). In the large- iV,- limit, however, 
all branches that have a gap in the dispersion relation at p = Q give contributions to P/T^ that are exponentially 
suppressed for large- iV^. Also for more general actions Eq. (14) thus is the starting point for the analysis of cut-ofi^ 
effects in the large-AT^ limit. In order to analyze power-like corrections to bulk thermodynamics in the continuum 
limit it suffices to analyze the properties oi E{p) = Ei{p). 

B. Rotational invariance, improved dispersion relations and improved thermodynamics 

For both, standard staggered and standard Wilson fermions the dispersion relation receives O(a^) corrections, 

3 

E\p) = / + 0{aV) + 0{a^Y.Pt) - (19) 

with p^ = p^ = X]fc=iPfc- Improving the dispersion relation moves the leading 0(0,^) lattice artefacts to some higher 
order n. For this purpose it suffices to construct an action with rotational invariance maintained to this order. To 



Wc give here relations for f^a = 0. The chemical potential can, however, be reintroduced at any stage through the substitution 

ap4 — > ap4 — ifJ-a. 
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sec how this works let us briefly discuss the dispersion relation for standard staggered fermions which is identical to 
Eq.(13). Expanding this relation for small values of apf^ and aE up to terms of order p^, E'^ yields 



(20) 



which has the solution 



E'{p)=p'-laV-lj:^'pt , (21) 

fc=i 

i.e. corrections to the continuum result, E{p) = p, start at ©(a^). In a similar way we may proceed for improved 
actions which have a Euclidean propagator that is rotationally invariant at 0{p*), i.e. it still may contain corrections 
at 0{a^), however terms proportional to are absent. The important new feature is that the low momentum 



V- /J " i'--! ^nr-fj, J- 

expansion of D, for the branch Ei surviving the continuum limit, now has the form 

(S? - /)[! + f{apk, aEi)] + OiaSD = 
with /{apk, aEi) = const a'^{Ef — p^) which leads to 

Ef=p' + 0{aV,) . 



(22) 



(23) 



In general, if the propagator can be factorized, up to order a"p^+^, into a factor E^ — p^, with Ei corresponding to 
the lowest root, times 1 plus some function f{apk,aEi) which starts at ©(a^) then 



a-'D{E,,p) = {Ef-p' + 0{aV^+')) [1 + f{apk,aE^)] + 0{a^pl+') 

and the dispersion relation receives C'(p"+^) corrections only i.e. it is ©(a""^) or 0{p") improved. 
Suppose now that the dispersion relation is 0{a"~^) improved, 



E-' 



p 



■0{a-pl+'). 



In this case the corrections to Ei/T start at 0{N^ "■), 



(24) 



(25) 



(26) 



where A is a polynomial with only even powers in 1 /N^ . 

We may now use this expansion to evaluate corrections to the continuum result for the press ure. To do s o we 
rearrange the argument of the logarithms in Eq. (14); with yi = Pi/T = N^api, y = p/T = -y/y^ + + 2/1 and 
A = A{yi/Nr, y^/Nr, ys/N^) one has 



1 + Ae-^i/^ = (1 + Ae-y) (l - S (l 



)) , 



(27) 



where B is given as B{y, A) — 1/{A ^e^ + 1) and A = z or z ^, respectively. After expanding in A one then finds 
for the logarithms appearing in Eq. (14) 

In (1 + yl exp {-NraE{p))) = In (l + Ae'^) - By A + 0{A^) . (28) 

While the first term will give the continuum ideal gas result for the pressure of a fermion gas, the second term and 
further higher order corrections will generate a systematic expansion in powers of l/Nr- This shows that improving 
the dispersion relation to a certain order in 1/Nj- immediately leads to an improvement of the high temperature limit 
of free energy, pressure and other thermodynamic observables derived from these to the same order. This also is the 
case for derivatives of P/T^ with respect to the chemical potential, such as the quark number density or higher order 
susceptibilities. We will expand on this in the next subsection. In the following it also will be convenient to introduce 
polar coordinates for y^, which allows us to think of the cut-ofl[ dependence entering only in the radial component, i.e. 
we take the arguments of A to be y/Nr, (p, 9 and rewrite it as 



A = A{y/Nr,cl>,e)=Y.^2k{<f>,e)ij-\ . (29) 

fc=i ^ 



The lowest non- vanishing contribution will appear at k = n/2. 
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C. Cut-off dependence of the pressure 



We are now in a position to discuss the systematic expansion of the pressure in inverse powers of Nt and analyze 
systematically cut-off effects and there dependence on temperature {aT = l/N^.) as well as chemical potential (/u/T = 
liaNT-). We insert Eq. (27) into Eq. (14) and obtain 

P f P 



y4 \ Ji4 



SB 



/ y^dy sinM6' / dcj) [\n (l - B{y, z) {l - e~y^)) + \n {l - B{y, z'^) {l - e-^^))] (30) 
Ja Jo Jo 



(27r)3 

where the leading order ideal gas term is given by 
P 



2^4 



1 f 

= — / dyy2 [In (l + ^ e"^) + In (l + ^"^e"^)] 
^ Jq 

47r^ n A ^ . A* \\ 77r2 30 / ^ 15 / a* V 



SB 

2 



(31) 



Note that the well-known ideal gas result for the pressure is an even Bernoulli polynomial of degree 4 in the 
chemical potential. In the following we will argue that similarly to the leading term also the corrections at order 
-A^~^" are Bernoulli polynomials in /i/(7rT) of degree (4 + 2n) which might be less obvious. 

As the entire cut-off dependence in Eq. (30) arises from the cut-off dependence of the dispersion relation, i.e. from 
A = A{y/Nr,4',9), we may first expand the pressure in a Taylor series around yA{0,(p,6) = 0. After noting that 
derivatives with respect to yA can be replaced by derivatives with respect to y we arrive at 



1 1 /" 

V^^VsB ~ 4^§^-/o "^'Jo" Jo ^"^'""^ V ar-^ ' 

Using the representation of A given in Eq. (29) we now can rearrange the above series and obtain an expansion in 
even powers oi 

^-{^)^^=tA^M^){i^y\ (33) 

where results from the integration over angular variables and the chemical potential dependent part, P2k {i^/t^T), 
is given in terms of a Bernoulli polynomial that results from the integration over the radial coordinate y, 

--Y.i-^r^-t\^r^^{^r'-'' (34) 



21 ;23+2'=-i64+2fe Vttt; 



with bn denoting Bernoulli numbers. The functions P„ have been normalized such that P„(0) = 1. Note that the ideal 
gas result can similarly be written as AqPq{pl/'kT) with Aq = 77r^/180. We give some more details on the derivation 
of Eq. (33) and the definition of ^2fe in Appendix A. The point we want to stress here is that the dependence of 
cut-off effects on the chemical potential arises only through a Bernoulli polynomial of degree 4 + 2k. For the particular 
complex argument appearing in Eq. (34) the Bernoulli polynomial contains only even powers of the chemical potential, 
reflecting the particle anti-particle symmetry, /i <-> —fi, of the partition functions. Wc also note that all terms in 
Eq. (34) are positive. Cut-off effects thus become larger with increasing /i/T, although the effect is small as exact 
evaluations of the pressure for fixed show [9] . 

Eqs. (33) and (34) are the main result of this section. With this we have shown that in the high temperature, 
ideal gas limit of QCD cut-off effects in the fermion sector of bulk thermodynamic observables can be traced back to 
the cut-off dependence of the dispersion relation for the quark propagator. The structure of the cut-off dependence 
is preserved at non-zero chemical potential fi, i.e. actions which are improved to a certain order in aT = 1/Nt at 
11 = are so also for /x > 0. Moreover, we find that the /x-dependence of cut-off effects is universal, i.e. /x-dependent 
correction factors at 0{N") are independent of the discretization and improvement schemes and are proportional to 
a Bernoulli polynomial of degree 4 + n. 

In the following we give explicit results for cut-off dependent corrections to the continuum ideal gas limit for several 
fermion actions and compare the asymptotic large N^. behavior derived in this section with an explicit evaluation of 
partition functions at small values of iV^ . 
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IV. STAGGERED TYPE QUARKS 

We now want to apply the general considerations presented in the previous section to the case of staggered fermions. 
Wc restrict ourselves to the class of staggered actions which contain terms with fermion and anti-fermion fields 
separated by up to three links, 

Sf = Y,x{x)M{x,y)x{y) , (35) 

x,y 

with 

I 

M{x,y) = ^V^iix) \ ^ c^^o[d{x + ifi,y)-S{x-ifi,y)] 

tJ.=l \i=l,3 

+ X! X! ^^'^ [6{x + fi+ jv, y)-5{x- jl+ ju, y)] \ . (36) 

To reproduce the correct continuum limit the expansion coefficients are constrained by the relation 

ci,o + 3c3,o + 6ci,2 = ^ . (37) 

Within this class of staggered fermion actions one can construct actions that are rotationally invariant up to 0{p'^), 
i.e. the denominator of the fermion propagator, D{p), becomes a function of only up to this order. This can be 
achieved with the additional constraint [3] 

ci,o + 27c3,o + 6ci,2 = 24ci,2 -^i:^^ ci,2 = 1/48 + C3,o • (38) 

Particularly well-known versions of 3-link actions that are currently used in studies of QCD thermodynamics with 

staggered fermions are the Naik action [2], for which ci^2 = 0, and the p4 action [3] where c^^ = 0. 

For the general three link action, not yet restricted to the subset of rotational invariant propagators, the propagator 
in momentum space is a polynomial of degree 3 in sin^(ap4), i.e. it is of the generic form given in Eq.(15). Following 
the discussion presented in the previous section one may work out the dispersion relation on the branch that survives 
in the continuum limit. For the leading correction to the continuum result one finds 

^ = 1 - i(l - 48(ci,2 - C3,o)) {p' + i^yyt]o:' + 0{a^) , (39) 

where wc have used Eq. (37) to eliminate ci o- 

As corrections to the continuum dispersion relation, Ei = p, start at ©(a^), bulk thermodynamic observables will 
receive cut-off dependent corrections starting at the same order, (aT)^ = l/N^- Eqs. (39) and (38), however, show 
that all three link actions with a rotationally invariant propagator will lead to an 0{a'^) improved dispersion relation 
and thus also to O(a^) improved bulk thermodynamics. In the class of three link staggered fermion actions one finds 
for the O(a^) improved dispersion relations. 




1= 

-^(1 + 1344c3,o)) - i Y.P'p')j + l^(^'iP2P3)^a 



+ 0(a«). (40) 



We note that all three link actions lead to identical leading order corrections to the continuum dispersion relation, 

but differ at subleading order. According to the discussion presented in the previous section, this feature will carry 
over also to the cut-off dependence of bulk thermodynamics at zero and non-zero chemical potential. 
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action 




A4/A0 




standard staggered 


248/147 


635/147 


3796/189 


Naik 





-1143/980 


-365/77 


p4 





-1143/980 


73/2079 


standard Wilson 


248/147 


635/147 


13351/8316 


hypcrcube 


-0.242381 


0.114366 


-0.0436614 


overlap/ 
domain wall 


248/147 


635/147 


3796/189 



TABLE I: Coefficients in tlie series expansion (Eq. (33)) for the pressure calculated using different actions. The expansion 
coefficients have been normalized to the ideal Fermi gas value at vanishing chemical potential, = {p/T'^)sB{n/T = 0) = 
77rVl80. 



Wc have evaluated the leading cut-off dependent corrections to the continuum ideal gas pressure for a generic 3-link, 
staggered fermion action up to 0{1/N^) using Mathematica to perform the expansfon of the integrand appearing in 
Eq. (30) in inverse powers of Nt- For the pressure we find, 



180 
7^ 



P 

rp4 



P 

1J14 



SB. 



248(1 -48(ci,2-C3,o)) f ^ 

147 ^ Wry I Nr 



^ /635 ^ 16256(c3,o - 01,2) ^ 2365248(c3,o - ci.z)^ \ f j^-^ f 



+ 



\U7 
73 
2079 



35 



245 



KnTJ \N. 



(1 + 6528c3,o + (1 - 48(ci,2 - C3,o))(571 + 1018368c?_2 



+ 51888c3,o + 13317120^ - 48ci,2(1217 + 48960c3,o))) Pe (^) (-^ 

+0((7r/iV.)8) , 



(41) 



with Pn defined in Eq. (34). Results for the standard staggered (1-link) action as well as for the class of C(a^) 
improved 3-link actions, for which ci,2 = 1/48 + 03,0, are easily read off from Eq. (41). In Table I we give the 
expansion coefficients for the standard staggered action as well as the Naik and p4 actions. Wc note that coefficients 
of the sub-leading corrections to the pressure calculated with the standard staggered action as well as with the Naik 
action are larger than the leading expansion coefficient. The subleading correction for the Naik action is a factor 137 
larger than for the p4-action. 

It is instructive to compare the asymptotic large Nr behavior of the cut-off dependent corrections to the continuum 
ideal Fermi gas with the exact results for the pressure evaluated at fixed values of N^- by using Eq. (6). At vanishing 
chemical potential this is shown for standard staggered, Naik and p4 improved fermions in Fig.l. Here we have 
normalized the exact lattice results to the continuum Stefan-Boltzmann limit, 77r^/180. We note that the exact 
calculation of the partition function at fixed vahics of reproduces the basic pattern derived from the leading orders 
of the large- expansion. Cut-off effects in the pressure are reduced when using improved discretization schemes. 



V. WILSON TYPE QUARKS 



We will now consider the case of Wilson fermions. Since in the interacting case cut-off effects arise in the Wilson 
fermion formulation already at 0{ag^) improvement schemes developed for Wilson fermions have primarily focused 
on removing terms linear in the cut-off (clover improved actions [4]). Actions that are 0(0^3°) improved received 
much less attention and, in fact, are currently not exploited in studies of QCD thermodynamics with Wilson fermions. 
Nonetheless, removing or at least reducing 0{a'^) effects in the Wilson fermion formulations is as important as in 
the staggered case when one wants to perform studies of QCD thermodynamics at high temperature within this 
discretization scheme. We will discuss here a generic Wilson type action with couplings constrained to a hypercube 
of size {2ay. This type of action is discussed in the context of truncated fixed point actions [6, 15]. Some exploratory 
studies of thermodynamics with dynamical hypercube fermions have been presented in [16]. The Dirac matrix for 
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0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.00 0.02 0.04 0.06 0.08 0.10 

FIG. 1: The pressure for staggered type quarks normalized to the continuum Stefan Boltzmann limit. The left hand figure 
shows results for the (standard) 1-link action and the left hand figure gives results for the Naik and p4 actions. In both cases 
we compare exact results, calculated on lattices with temporal extent Nt and infinite spatial extent, with the leading 1/N" 
correction. Note the different vertical scales. 

this generic Wilson action can be written as 

4 

M{x,y) = J2^^pA^~y) + Kx-y) (42) 

with 

Pf,{x-y) = pi [S{y,x + fi) - S{y,x - fi)] 

+P2^[S{y,x + fi + D) - S{y,x ~ ii + i))] 

+P'i^[6{y, x + P + + p)- 6{y, x~ + + p)] 

+P4 ^ [5{y, X + p. + + p + a) — 8{y,x — (i + V + p + a)] (43) 

for the vector terms and 

K^-y) = ^0 5{y,x) 

+Ai ^[(5(y, x + fj:)+5{y,x- fi)] 

+ A2 ^ [(5(j/, .T + /i + i>) + (5(y , X - /i + z>)] 

+^3 X! ['^(^'^ + A + + p) + S{y,x- P + V + p)] 

+A4 ^ [5{y,x + ji + v + p + a)+6{y,x- il + v + p + cr)] (44) 

for the scalar ones. The sums over v, p, a extend over positive and negative directions and are mutually orthogonal to 
each other and to p. Examples for this type of action are the standard Wilson action (including the clover improved 
version of it) and the hypercube truncated perfect action [6, 15], with coefficients as listed in Table II for the massless 
case. 

The denominator of the fermion propagator, D{p), is most conveniently written as 

D{p,pi) = R{p) + 2P{p) cos{api) + Q{p) coa^{api) . (45) 

Note that this complies with Eq.(16) since cos(ap4) is to be rewritten as 1 — 2sm^{api/2). The coefficients P, Q and 
R depend on p and are listed in Appendix B. The requirement that in the limit of ap^; for all p — 1, ..,4 the 
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Hypercube 


Wilson 




Hypercube 


Wilson 








A„ 


1.852720547 


4 


pi 


0.136846794 


1/2 


Ai 


-0.060757866 


-1/2 


P2 


0.032077284 





A2 


-0.030036032 





P3 


0.011058131 





A3 


-0.015967620 





P4 


0.004748991 





A4 


-0.008426812 






TABLE II: Coefficients pi and Ai for standard Wilson quarks and for the hypercube action. 



continuum dispersion relation is obtained leads to the constraints (see also [17]) 

2pi + 12p2 + 24p3 + 16/34 = 1 , (46a) 
Ao + 8Ai + 24A2 + 32A3 + I6A4 = . (46b) 

Moreover, demanding that the so-called doublers receive a non-vanishing mass at the corners of the Brillouin zone, 
= IT for one or more component /U, introduces further conditions to be satisfied by the A coefficients, 

Ao + 4Ai - I6A3 - I6A4 (47a) 

Ao - 8A2 -I- I6A4 ^ (47b) 

Ao - 4Ai -I- I6A3 - I6A4 ^ (47c) 

Ao - 8A1 -I- 24A2 - 32A3 -I- I6A4 ^ (47d) 

In the standard Wilson action these conditions are realized such that the so-called doubler masses arc 2/a, 4/a, 6/a 
and 8/ a while for hypercube fermions these masses are ~ 2/a with deviations of at most 2 per cent for aU corners ^. 

We may now go through the steps outhned in Section II to determine the cut-off dependence of dispersion relations 
and bulk thermodynamics. Wc start by expanding the quark propagator, Eq. (45), for small momenta, taking into 
account the constraints given in Eq. (46), 

D{p,Pi) = {paf - i(-l + 12p2 + 48p3 + 48p4) ({paf - Y^'^P^)^ + ^(«') ' (^8) 

With the couplings in the hypercube action as listed in Table II the propagator receives 0{a^) corrections, which 
however are small when compared to the standard 1-link Wilson action. For the latter the prefactor in Eq. (48) equals 
1/3 while for the hypercube action it is —0.0479. Of course, with an appropriate choice of the couplings p2, p3 and 
Pi one can also eliminate the 0{a?) corrections to the propagator of the hypercube Wilson action completely, making 
it rotational invariant up to that order. 

Turning now to the dispersion relation, from Eq.(45) one obtains, for the branch relevant to the continuum limit. 




sinh^^::|ij=--|^i+ - ' .^^ ) . (49) 

The low momentum expansion of this dispersion relation yields, 

^ = 1 + + 12p2 + 48p3 + 48p4) ( / + \ ^/fc ) + 0{a^) . (50) 

Wc note that for standard Wilson fermions (p2 = /O3 = /54 = 0) the leading order correction to the continuum 
dispersion relation is identical to that of staggered fermions. This will, of course, carry over also to the leading 
corrections to bulk thermodynamic observablcs. Moreover, a choice — 1 + 12/92 + 48/03 + 48p4 = not only removes 
the rotational symmetry breaking term, Y^k=i'Pk-> propagator and the dispersion relation, but eliminates the 

0{a?) corrections completely. 



^ Note that with five equations the A coefficients can be determined uniquely once the doubler masses have been fixed. Indeed, fixing 
these masses to 2/a leads to the HF scalar couplings within 2 per cent accuracy. 
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FIG. 2: Pressure for standard Wilson fermions (left) and liypercube Wilson fermions (right) normalized to the continuum 
Stefan Boltzmann result. In the left hand part we compare the results obtained with standard Wilson fermions to that of 
standard staggered fermions. Both discretization schemes have identical cut-off errors in leading order l/N'i. The right hand 
part compares results for the hypercube action with the p4 action. 



Using the dispersion relation for the hypercube action, Eq. (49), in our basic relation for the pressure, Eq. (30), we 
again can perform a systematic expansion for large N^-. With this we obtain for the leading cut-off corrections to the 
continuum pressure. 



180 
7^ 



P 



P 



SB. 



48p3 + 48p4)P2 



127 
'735 



25 + Mmpl -f 55872p^ - 2148p4 + 55872p^ 



+2328p3(-l -i-'iO/"4 
\ttTJ 



3p2(-209 + 9312p3 + 9312p4)) Pi (^) 



TT 



n 



+ O((^/iV0«) 



(51) 



The explicit form of /ig is rather lengthy. It is given in Appendix C. 

In Fig. 2 we compare exact results for the pressure obtained with standard Wilson fermions and the hypercube 
action with the leading order large- A't- expansion obtained in this section. As can be seen in Fig. 2(left) the standard 
Wilson action and standard staggered fermions have identical leading order cut-off errors and consequently follow an 
identical pattern for large Nr, e.g. for A'r^lO. Cut-off effects for the hypercube action are significantly reduced and 
of similar magnitude as for the improved staggered actions (Fig. 2(right)). Cut-off dependent corrections at small 
values of Nr are already close to the result obtained from the leading order expansion reflecting that higher order 
corrections are small. This also is apparent from the higher order corrections listed in Table I. 



VI. OVERLAP AND DOMAIN WALL FERMIONS 



Let us now comment on bulk thermodynamics with chirally invariant fermion formulations, e.g. overlap [18] and 
domain wall fermions [19, 20]. The free field limit of chiral fermion formulations [21] has recently been used to 
analyze the cut-off dependence of spectral functions at high temperature [22] . We follow this discussion here closely 
to explore the cut-off dependence of bulk thermodynamics in lattice formulations that use chiral fermion formulations, 
e.g. overlap or domain wall fermions with standard Wilson kernel. 

It is known that in the free field limit the dispersion relations for overlap and domain wall fermions are given by 
the standard (staggered) dispersion relation up to a constraint, which only alters the large momentum part of the 
dispersion relations [21] . As we have seen in the previous sections this part of the spectrum does not influence the large- 
N-r expansion of thermodynamic quantities. Following the general discussion of cut-off effects given in Section II we 
thus expect that the 1/Nr expansion for the cut-off dependence of the pressure and related thermodynamic quantities 
is identical to that of standard staggered fermions. As outlined in Section II this also includes the dependence on a 
non-zero chemical potential. Nonetheless, in particular for the case of overlap fermions the question has been raised 
whether the fact that one has to take a square root of e.g. the Wilson fermion operator may lead to difficulties in 



12 



reproducing the correct continuum thermodynamics at non-zero chemical potential [12]. Indeed, the presence of a 
square root in the basic relation for the partition function (see Eq. (57)) leads to cuts which so far did not show 
up in our general discussion of bulk thermodynamics in the ideal gas limit. We will, however, show that the cut 
contributions are exponentially suppressed and do not contribute to the 1/Nt expansion of the pressure. 

We outline here the basic steps that allow to relate also in the overlap formalism the basic properties of the dispersion 
relation with the cut-off dependence of the pressure. Some more technical considerations are given in Appendix D. 

We work with the following form of the (massless) overlap operator [21, 22], 

D = -(l + ^=\. (52) 

Here X = Dw — 1, where Dy/ is a nonchiral, undoubled fermion operator like the standard Wilson operator discussed 
in the previous section. In momentum space X{p,pi) takes on the form 

X{p,Pi) ==- [i^lu sin(ap^) + 2 ^ svo?{ap^/2) - 1 J • (53) 

\ u—l v—1 / 

In order to keep the notation simple we do not explicitly introduce a chemical potential. It can, however, at any step 
be re-introduced through the replacement api — > api — ifxa. 

The fermion propagator of the overlap operator is then given by 

D{p,p,) ='-(l- S^illd^l , (54) 

where 

4 4 

A=A{p,p4) = l-2Y^ sin^ (ap./2) = {p, P4) =Y1 sin' {ap.) • (55) 
Prom Eq. (54) we see that the poles of the fermion propagator are given by 



V^42 + 52 -A = ^ = &n<l A>Q . (56) 

Thus, overlap fermions have indeed exactly the same dispersion relation as naive fermions upto a constraint, ^ > 0. 
The constraint is automatically satisfied at small momenta but has no solutions near the edge of the Brillouin zone, 
thus avoiding the doubling problem. As pointed out in [22] this carries over to similar relations for the propagator 
and the poles of free domain wall fermions in the limit of infinite 5th dimension. 
With this we obtain for the pressure 



Pa* = 



r d^ap-^ y In f 1 - , } i . (57) 

UP4 



In Eq.(57) the appearance of a square root, of course, does not create any problems; its argument is positive and the 
pressure thus can be evaluated for any fixed temporal lattice extent A'',-. The discussion of the large- A^ expansion, 
however, becomes somewhat more complicated than in the cases considered in the previous sections. We also note 
that we did not reduce the integration range, as it has been done for naive fermions to eliminate contributions from 
doublers. Integrals and momentum sums still extend over the whole Brillouin zone. Let us denote x = sin(ap4/2). 
The expression for the pressure may be shown to take the form 



Pa''=—— d-^ap— > In 1 — , ' \ , (58) 

(27r)3i_^ ^Nr^ \ V2 Vc2x2 + 62 ' ' ^ > 



where 

I? = y s\n^ api/2s\n^ apj/2 + l/S,= \{l+w'^ +A{c^ - c^)) , 



l<0 
3 



C2 



= y^sin^ api/2 , uP' = sin^ apj . (59) 
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The argument of the logarithm may be written as a quadratic form in ^ = V c^x^ + 6^, given by 
(C^ + 0^72^ + {c^{c^ - 1/2) - fe^)) /(V2c^0- Hence the numerator may be factorized to obtain 



Pa^ = 



j d^ap^ (ill (Vc^x'^ + 62 - ^+ j + In (^y/c^x^ + 62 



, , (60) 

with ^± = -c2/\/2± v/62 - (c4 - c2)/2 = -cVv^± i/(l +a;2)/8. The last term in Eq. (60) is easily identified as 
a vacuum term that gets canceled once we subtract the zero temperature contribution to the pressure. The third 
term has a structure similar to that discussed for staggered fermions. It corresponds to a dispersion relation that 
approaches a non-zero value for small momenta and thus will not contribute to a 1/Nr expansion for the pressure, 
i.e. it only leads to exponentially small contributions. We thus have to understand the Nr dependence of the first 
two terms. To do so we first carry out the sum over p4 which can be done in close analogy to the cases discussed 
above by first considering the sums 

^ ap4 ap4 * ^ 

This sum can be evaluated using the same contour integral technique [14] exploited in the case of staggered and 
Wilson fermions. It, however, needs a bit more care as through the replacement of = sin2(ap4/2) by a complex 
variable, = —{z — 1/zY /A, the square root in Eq. (61) will have cuts. These contribute to the contour integral. 
Additional residues appear when doing tlic^ sum for ^+ as poles appear in Eq. (61) for small momenta when ^+ > 0. 
We discuss the calculation of S"^ in Appendix D and show there that contributions arising from the cuts differ from a 
vacuum contribution only by exponentially small terms and thus do not contribute to the large- A^t- expansion of the 
pressure. This expansion thus receives contributions only from the first term in Eq. (61) and only from that part of 
the momentum integral where > 0. We thus find, 

^ ^ [Pa4 _ (p^4)^] ^4 = 1^ 1^ A^ap J2 111 (Vc2sin2(ap4/2) + 62 - e+) 



2Ari 



(27r)4 



^ d^ap In (^c2sin2(ap4/2) + 62 - + 0(e-"^- ) 



(62) 



where d{p) denotes the section of the 3-d Brillouin zone, in which ^+ > 0, i.e. for which 1 -|- a;2 > 4c'' holds. As 
discussed in Appendix D this result also holds for non-zero chemical potential, with the replacement ap4 — > ap4 — i^,a. 

We are thus left with an expression very similar to what we have discussed in the previous sections and can express 
the pressure for overlap fermions through the poles of the denominator appearing in Eq. (61), 

sinh2(£+/2) = 1 (z+ - ^) ' = ^^1+^) ^ ^iiili' E+=uj^ . (63) 

We now can integrate Eq. (61) to arrive at explicit expressions for Matsubara sums of logarithms. After subtracting 
again the zero temperature contribution to the pressure one then obtains for the pressure in the overlap formalism 
an expression that agrees with the corresponding result for staggered fermions, Eq. (14) in a limited region of the 
Brillouin zone that includes p = 0. Aside from this the pressure for overlap fermions receives only contributions that 
are exponentially small in N^.. This confirms that the large- A''^- expansion for the pressure in the overlap formulation 
agrees with that of staggered fermions. As indicated at the beginning of this section this could have been anticipated, 
for we know that the two actions have the same dispersion relation at small momenta. 

For fixed values of the exact expression for the pressure, Eq. (57), of course, differs from that of naive staggered 
fermions and cut-off effects thus will be different. Using Eq. (57) we have evaluated the pressure for small values of 
Nt explicitly. The results are compared to corresponding calculations performed for standard staggered and Wilson 
fermions as well as to the common leading order large-A^^ expansion result in Fig. 3(left). We note that the cut- 
off dependence introduced by the overlap operator for small values of Nt is smaller than in the staggered fermion 
discretization scheme. In fact, it is remarkable that already for Nr = 8 the exact result for overlap fermions is close 
to the leading order result of the large- A^^ expansion. For A't->12 all 1-link discretization schemes yield very similar 
results and are in good agreement with the leading order large- A^^ expansion. 

We finally comment on the ideal gas limit for domain wall fermions. This has previously been analyzed numerically 
on finite space-time lattices [23]. In the analysis of the cutoff dependence of bulk thermodynamics with domain 
wall fermions an additional complication arises. One has to introduce additional degrees of freedom, Pauli-Villar 
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FIG. 3: Pressure for overlap fermions based on the application of the Ginsparg- Wilson operator on a standard Wilson operator 
(left) and for the domain wall operator using the same kernel with domain wall height M5 — l(right). The large- A^'t expansion 
for both discretization schemes coincides with that of the standard staggered action. In both figures we compare with the 
leading order result. 



regulators, which eliminate divergences that arise from the infinite number of heavy fermion fields present in the limit 
of infinite 5th dimension [24]. The calculation of the free fermion partition function for domain wall fermions, i.e. 
the free fermion determinant, thus is somewhat more involved [25]. Nonetheless, after diagonalization in momentum 
space the regularized domain wall fermion determinant takes on a simple form in the limit of infinite extent of the 
5th dimension. For massless fermions one obtains, 

detDow = YlJ^{p,P4) , (64) 

P,P4 

with 

M{p,Pi) -21+ , . (65) 

y ^(1 + ^2 + 14/2)2 _ 4^2 J 

Here has been defined in Eq. (55) and W = 1 — M5 — "lY^v^i siii^(aPi//2) for domain wall height Afs. 

As mentioned before also domain wall fermions have the same dispersion relation as staggered fermions, up to 
a constraint, which cuts off the spectrum for large momenta. The large- TV,- expansion of bulk thermodynamics in 
the ideal gas limit thus again will be identical to that of naive staggered fermions. Moreover, for vanishing lattice 
spacing in the 5th dimension and infinite extent in this direction the domain wall Dirac operator is identical to 
the overlap operator. It thus may not be surprising that the cut-off effects in the domain wall and in the overlap 
formalism are quite similar in the free field limit. Differences arise because the sign function is introduced in the 
domain wall formalism through a discrete set of additional fermion fields [24] . We limit our discussion of cutoff effects 
in thermodynamics with domain wall fermions here to a comparison of the exact evaluation of the pressure in the 
domain wall and overlap discretization schemes. With Eqs. (64) and (65) the pressure of free massless domain wall 
fermions is obtained as 

— EE [Pa^ - {Pa\] ^ j^J d^apY,\nU{p,p,) -J^J_ d^ap\n//{p,p,) . (66) 

ap4 

In Fig. 3(right) we show results from a direct evaluation of this expression for some values of Nr and for domain wall 
height A/5 = 1, which corresponds to the case discussed above for overlap fermions. As can be seen cutoff effects differ 
for small values of Nr but rapidly become similar for ~ 8. 



VII. CONCLUSIONS 



We have analyzed here cut-off effects that arise in lattice calculations of bulk thermodynamic observables at high 
temperature and vanishing as well as non-vanishing chemical potential. We have shown that in the high temperature, 
ideal gas limit of QCD cut-off effects in fermionic observables can be traced back to the cut-off dependence of the 
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dispersion relation for the quark propagator. The structure of the cut-off dependence is preserved at non-zero chemical 
potential /i, i.e. actions which are improved to a certain order in aT = l/N^ at = arc so also for /i > 0. We 
could show that the /i-dependence of cut-off effects is universal, i.e. /x-dependent correction factors at 0{N") arc 
independent of the discretization and improvement schemes. They are given in terms of a Bernoulli polynomial of 
degree 4 + n; cut-off corrections at order 1/N" depend on an even polynomial in /x/ttT of degree i + n with positive 
coefficients only. 

We have shown that the considerations on cut-off effects in thermodynamics with staggered and Wilson fermions 
carry over to chiral fermion formulations. So far the latter have been used mainly with standard Wilson fermion 
kernels. The discussion presented in this paper suggest that studies of bulk thermodynamics performed with chiral 
fermions could profit a lot from the application of improvements schemes similar to those used in studies of bulk 
thermodynamics with staggered fermions. This certainly can be achieved using the concepts that have been developed 
to improve standard staggered and Wilson fermion action [26]. 
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APPENDIX A: LARGE iV^ EXPANSION OF THE PRESSURE 

We give here some further details on the derivation of the series expansion of the pressure, Eq (33), and give explicit 
formulas for the expansion coefficients An as well as the polynomials P„ introduced in Eq.(33). 

In order to derive Eq. (32) from Eq. (30) as the result of a Taylor expansion yA it is helpful to note that 

dln[l-B{y,z){l-e-y'')] ^ -1 _ 

dyA z exp(j/) cxp(yA) -I- 1 

It then is obvious that taking further partial derivatives of this logarithm with respect to yA at yA = is equivalent 
to taking derivatives with respect to y at yA = 0. This yields Eq. (32). 

The series representation of the pressure, Eq (33), is then obtained from Eq. (32) after noting that the expansion 
coefficient at order l/N"^^ receives contributions from all terms in Eq. (32) with n < A: as A is a polynomial in 
y/N^ starting at 0{1/N^). Aside from a complicated factor that only depends on the angular variables {(p, 9), these 
contributions arc proportional to y^*^. Irrespective of the structure of the dispersion relation of a particular action 
under consideration the contribution of the radial integral in Eq. (32) to an expansion coefficient at order l/N"^^ thus 
has the unique form 

The angular contribution is given by 

L fc -1 /'TT n2TT J>- 

^^^n\jQ Jo (k,,..K)o=i 

where the inner sum extends over all n-tuples of positive even integers with — 2^- Here we have introduced a 

multiplicative normalization factor, the Bernoulli number 64+2*;; that arises from the integration of the radial variable 
y at vanishing chemical potential and normalizes the polynomials P2k{x) to unity for x = 0. 

For a few small values of the expansion parameter k we give here the explicit form of the polynamials P2k{x) defined 
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in Eq. (34): 

/ \ - 30 o 15/1 



7-^ ' 7 

n / N 1 147 2 10-5 4 21 g 
31 31 31 

620^2 _L 490^4 _L I'^'^^e _L 1^ a-s 



„ , , ^ 620 2 490 4 140 6 15 

12573 2 2046 , 462 ^ { 

X H X H X H a; i 

2555 511 365 511 255 



, , 12573 2 2046 , 462 . 99 o 33 /«.^ 
Peix) = iH H x^ H H H x^° (A4) 



APPENDIX B: THE PROPAGATOR FOR HYPERCUBE WILSON FERMIONS 

using the shorthand notation, c, = cos api and Sj = sin api we write the the denominator of the fermion propagator 
for Wilson type fermions, Eq. (45), as 

D{P4,P) = (ixr/C? +(5^) + 2cosA;4(iTr/Ci/C2 + K1K2) +cos2 A;4(iTr/Ci + -J2) (Bl) 

where 

Ki{p) = Ao + 2Ai(ci + C2 + C3) + 4A2(ciC2 + C2C3 + C3C1) + 8A3C1C2C3 (B2a) 

K2{p) = 2Ai + 4A2(ci + C2 + C3) + 8A3(ciC2 + C2C3 + C3C1) + I6A4C1C2C3 (B2b) 

5{p) = 2pi +4p2(ci + C2 + C3) + 8p3(ciC2 + C2C3 + C3C1) + I6P4C1C2C3 (B2c) 

and the /Cj matrices are given by 

3 3 

^1 = X] '^i'^i ^2 = ^ 7i A • (B3) 

i=l 1=1 

The coefficients appearing in this matrix representation are, 

ai{p} = 2si(pi + 2p2(c2 + C3) + 4P3C2C3) (B4a) 

a2{p) = 2s2(pi+2p2(c3+Ci)+4p3C3Ci) (B4b) 

a3{p) = 2s3(pi + 2p2(ci + C2) + 4p3CiC2) (B4c) 

Pl{p} = 4si(p2 + 2p3(C2+C3)+4p4C2C3) (B4d) 

(3^(0) = 4S2(P2 +2p3(c3 +ci) +4p4C3Ci) (B4e) 

Mp) = 4s3(p2 + 2p3(ci +C2) +4p4CiC2) (B4f) 

Finally, the expressions P, Q and R used in Eq.(45) and (49) arc obtained as 

P{p) = ^TrK.iK.2 + K1K2 Q{p) = ^Tr/C^ + 4-5'^ R{p) = ^Tr/C? + kI + 6^ (B5) 



APPENDIX C: THE EXPANSION COEFFICIENT he 

We here give the explicit result for the sixth order expansion coefficient he appearing in Eq. (51), which has been 
obtained using Mathematica. It reads 

—73 

he = —— (-2288 + 7344OOA3 + 47OOI6AI + (1175040 -2585088r+ 1410048/92)A3A4 
8316 

+82584p2 - 918432p^ + 288000p3 - 6846336p2P3 - 12690432p^ + 247392p4 
-6345216p2P4 - 23376384p3p4 - 10685952p^ + 29376A2(2 - Srf + 1836Ai(l - 2rf 

-1762560A^r - (940032 - 470016r)A^r + 1057536A3r^ + 5967r3 
+58752A2(-2 + 3r)(4A4(-l + r) + A3(-5 + 6r)) 

+14688Ai(-l + 2r)(4A4(-l + r) + A2(-2 + 3r) + A3(-5 + 6r))) (CI) 
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z - plane 




FIG. D.l: The locations of the poles and branch cuts of H{z) and the choice of contour for A^t ~ 4. The four poles of f{z) only 
exist for q > 0. Also shown is the unit circle (dotted line) and the location of poles arising from the exponential integration 
kernel introduced in Eq. (D3). 



where we have introduced the short hand notation, 

r = 8p2 + 32p3 + 32p4 



(C2) 



APPENDIX D: MATSUBARA SUMS FOR OVERLAP FERMIONS 



We give here some further details on the evaluation of Matsubara sums appearing in calculations with overlap 
fermions. Throughout this appendix we make the dependence on the quark chemical potential explicit. We want to 
evaluate the sum 



S{a) ^ ^ In (Vc2x2 + 62 _ 



(Dl) 



where x = sin((ap4 — i/ia)/2), and &, c and a are functions of the three-momenta {pi,P2,P3) and ap4 = 27r(n + l/2)/A^^ 
with n = 0, ±1, ±iNr - l)/2, Nr/2. 

To evaluate it, we first differentiate with respect to a. 



dS _ 1 ^ 
da ~ ~iv7^ 



1 



Nr ^ y/c'^x'^ + 62 - a 
Pi " 

This sum can be thought of as being part of the result of a contour integral over the function. 



H{z) 



1 



y^_c2 1)2/4 + 52 _ ^ ^ (^'~- + e^"^^- ) ' ziz^""- + e^-~^ ) 



(D2) 



(D3) 



where the second equality defines the function f{z). The function H{z) has poles and cuts; the sum given in Eq. (D2) 
is related to the poles on the unit circle which arise from the last factor in the denominator of Eq. (D3). We want 
to evaluate the contour integral along the path shown in Fig. D.l. As long as a = a{p) is positive there are four 
additional poles arising from f(z) which are located at ±zq and ±1/zq with, 



zo = e 



2 

and coshx = 1 + — (6^ - a^) 



(D4) 



The residues at these points are listed in Table D.l. 
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Pole: z 


Residue of H{z) 


Zo 




2a 


1 






- (1/^0)2 




+ 1 


-Zo 




2q 


1 




Z{) 


— ("1 /7n"l2 
U/ 20J 


\Zo^ ) ^ 


-1- 1 


1/20 




-2a 


1 




- 


- (1/20)2 




+ 1 


-1/zo 




-2o 


i 




. 

^0 


- (1/^0)2 




+ 1 



TABLE D.l: The four residues of ll(z) at the poles of ^(z) which exist for a > 0. 



Moreover, there arc branch cuts arising from the square root, which we take to be located in the intervals C~ = 
(—00, Wo), = ( — l/wo, 1/wo), C"*" = («;o, oo)j with sinh(wo) = hjc. We note that the contributions from cuts with 
Re2: > are identical to those with Rcz < 0. Moreover, one can map the contribution from the positive part of the 
central cut, (0, 1/wo), onto the interval C+. With this we arrive at the result, 



d5 ^ Aa f 1 1 

da c2 sinhx \e-''^T + 1 e^-(*-'''^) + 1 



■— Im / 

7^ Jw, 



Integrating this again we obtain 



(D5) 



S{a) = e(Q)-^ (ln(l + e-^-(^-'^'')) + ln(l + e-^-(*+'^'»)) + iV^x) 

+ -I-y^^ ( (,2e-.«);v. + i - (,-2e-.a)iv. + 1 ) + const. (D6) 

It easily is seen that the absolute value of the integral appearing here is bounded by const • exp(— iV^wo). When we 
insert S{a) in Eq. (62) we thus find that only the first term in Eq. (D6) contributes to the large- AT^ expansion of the 
pressure. 
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